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We study the intramolecular vibrational energy redistribution (IVR) dynamics of an effective 
spectroscopic Hamiltonian describing the four coupled high frequency modes of CDBrClF. The IVR 
dynamics ensuing from nearly isoenergetic zeroth-order states, an edge (overtone) and an interior 
(combination) state, is studied from a state space diffusion perspective. A wavelet based time- 
frequency analysis reveals an inhomogeneous phase space due to the trapping of classical trajectories. 
Consequently the interior state has a smaller effective IVR dimension as compared to the edge state. 



Investigating the dynamics of an initially localized vi- 
brational excitation of a molecule in terms of timescales, 
final destinations and competing pathways has been of 
considerable interest to chemical physicists for a num- 
ber of decades 01 0, [H, EL H, IS 0] • Due to the sustained 



theoretical [3, H J3 j3] and experimental efforts [3, H, 0, 0] 
it is only now that a fairly detailed picture of the 
intramolecular vibrational energy flow is beginning to 
emerge. Recent studies (il. fsl. Hol. suggest that IVR can 
be described as a diffusion in the zeroth-order quantum 
number space (also known as the state space) mediated 
by the anharmonic resonances coupling the zeroth-order 
states. The state space approachJ6j makes several predic- 
tions on the observables associated with IVR. Foremost 
among them is that an initial zeroth-order bright state 
| b) diffuses anisotropically on a manifold of dimension D 
much smaller than N (or N — 1 with energy constraint). 
N is the number of vibrational modes in the molecule. 
As a result the survival probability Pb(t) exhibits power 
law behaviour on intermediate time scales 
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with tif, = 5Z a |(&|a:)| 4 DCm g the dilution factor of the 
zeroth-order state \b) and \a) denoting the eigenstates of 
the system. Wong and Gruebele[I2] explained the power 
law behaviour from the state space perspective by pro- 
viding a perturbative estimate for D as: 
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with n — \1c\4 — n&| being the distance, in state space, 
from the state \b) to other states \i) and the sum is over 
all states \i) such that | — n& | < n. The zeroth-order 
quantum numbers nk are associated with the state \k) 
and the symbol A indicates a finite difference evaluation 
of the dimension due to the discrete nature of the state 
space. In practice one chooses two different distances n 
in the state space and evaluates eq.[5]and thus D ss D(n). 
The quantity 
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is a measure of the number of states locally coupled to 
| b). The difference in the zeroth-order energies is denoted 
by AE? b = E° - E® and V ib = (i\V\b). Notice that in the 
strong coupling limit, AEf b -C Vib, La, ~ 1 whereas in the 
opposite limit Lib ~ 0. Thus D v (n) can range between 
the full state space dimension and zero 12| . For further 
discussions on the origin and approximations inherent to 
eq. [5] we refer the reader to the original reference [12j. In 
the context of the present study it is sufficient to note 
that D oc Ni oc which has been confirmed in the earlier 
work[I2. 

Clearly Eq. [5] explicitly includes the various anhar- 
monic resonances and hence not only the local nature 
but the directionality of the energy flow is also taken 
into account. The main point is that the above estimate 
for D, which can be obtained without computing the ac- 
tual dynamics, is crucially dependent on the nature of the 
IVR diffusion in the state space. However, to the best of 
our knowledge, precious little is known about the dynam- 
ics associated with the state space diffusion. Our moti- 
vation for investigating the IVR dynamics in the state 
space has to do with the observation that the state space 
model shares many of the important features found in the 
classical-quantum correspondence studies of IVR[l|, 0,0). 
Classical dynamical studies identify the nonlinear reso- 
nance network as the crucial object. On such a network, 
directionality and the local nature of IVR arises rather 
naturally mainly due to the reason that molecular phase 
spaces are mixed regular-chaotic even at fairly high ener- 
gies. How does the mixed phase space influence the IVR 
diffusion in the state space? Is there any relation, and 
hence correlation, between the classical resonance net- 
work and the IVR dimension D in the state space? Is it 
possible that local dynamical traps in the classical phase 
space can affect the validity of Eq. [^]? Answers to these 
questions can have significant impact on our ability to 
control IVR and hence reaction dynamics. The issues in- 
volved are subtle and this preliminary work attempts to 
address the questions by studying a specific system. 

Although detailed classical-quantum correspondence 
studies of IVR have been performed [J, 0, 0, [3| on sys- 
tems with two degrees of freedom, in order to address 
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the questions one needs to analyze atleast a three degree 
of freedom case. This is due to the fact that the scal- 
ing theory of Schofield and Wolynes posits D = 2 as the 
critical scaling i.e., near the IVR threshold^. Thus for 
systems with two degrees of freedom the separation of 
diffusive and critical regimes is not very sharp How- 
ever, studying IVR from the phase space perspective is 
difficult in systems with three or more degrees of free- 
dom. In this study we use a time-frequency technique 
proposed by Arevalo and Wiggins [3] to construct a use- 
ful phase space representation of the resonance network 
for three degrees of freedom. Such an _approach, as seen 
below and in many recent studies TH, 16| , is particularly 
well suited for our purpose. Thus we choose an effec- 
tive spectroscopic Hamiltonian [I?) describing the energy 
flow dynamics between the four high frequency modes 
of CDBrClF. The Hamiltonian H = H Q + V res with the 
anharmonic zeroth-order part 
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has various anharmonic resonances coupling the four nor- 
mal modes denoted by s (CD-stretch), / (CF-stretch) 
and (a, b) (CD-bending modes) 
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The harmonic creation and destruction operators for the 
mode are denoted by a- and ou respectively. Note 
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that despite having four coupled modes the system has 
effectively three degrees of freedom due to the existence 
of a conserved polyad N = v s + (vf + v a + %)/2. In 
this work we choose N — 5 for illustrating the main idea. 
Similar results are seen in other systems and the details 
will be published later. The values of the various param- 
eters are taken from the fit in reference 17J (fourth col- 
umn, Table VIII). The Fermi resonance strengths k S j m 
are larger then the mean energy level spacings (13.7 cm -1 
~ 2.5 ps) of H for N — 5. Thus this is an exam- 
ple of a strongly coupled system and the multiple Fermi 
resonances render the classical dynamics irregular. We 
investigate the IVR dynamics out of two nearly isoen- 
ergetic zeroth-order states \v s ,Vf,v a ,Vb) = |5000), and 
|3301) denoted for convenience as |1) and |2) respec- 
tively. The experimentally accessible state |1) has en- 
ergy E° w 10571 cm -1 whereas the combination state 
has E° ~ 10567 cm -1 E°. We restrict our study to 
these two states although there are other close by states 
within a mean level spacing. In terms of their location in 
state space the state |1) is an example of an edge state 
whereas |2) is an example of an interior state. 
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FIG. 1: Survival probabilities Pb(t) and time-smoothed sur- 
vival probabilities Cb(T) for the two isoenergetic states. The 
left and right top panels show Pb(t) for |1) and |2) respec- 
tively. Intermediate time power law behaviour is exhibited 
by Cb(T) in both cases, dashed line (state |1)) and solid line 
(state |2>). Note that the C b (T) data is shown for T = 10 
ps indicating small oscillations about the respective dilution 
factors Ob- 



Given the Hamiltonian and the resonant couplings the 
number of states coupled locally to |1) and |2) can be 
estimated as Ni oc = 1.8, and 3.0 respectively. Combined 
with the fact that N e ff = a^ 1 sa 12 one expects fast 
IVR from both the states. Since N loc (\l)) < N loc (\2)) 
the decay is much faster at short times for |2). This 
is confirmed in Fig. [1] which shows Pb(t) for the states. 
However Fig. [T] also shows the time-smoothed survival 
probability [18J, [l9( 
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associated with the states and, importantly, highlights 
a power law behavior of Cb(T) ~ j^-d/i a t intermediate 
times - a sign of incomplete IVR[7, 12] . Note that the per- 
sistent recurrences in P\ (t) occur for much longer times as 
evident from the results for Cb (T) . Earlier works 18 , 1|| , 
in an apparently different context, have associated the 
power law behaviour with the multifractality of the eigen- 
states and the local density of states. The power law ex- 
ponent or effective IVR dimensionality in the state space 
are determined to be d\ ~ 1.8 and cfo ~ 1-3 which are 
smaller than the three dimensional state space. Inciden- 
tally, a purely exponential decay of P b {t) would imply 
Cb(T) ~ T^ 1 irrespective of the dimensionality of the 
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state space. More surprising observation from Fig. Q] is 
that the interior state shows faster short time IVR but 



(2) 

at longer times, despite 
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explores an IVR 



manifold of smaller dimension as compared to the edge 
state. Infact based on Ni oc and the strong couplings one 
would infer the opposite from Eq. [21 

Since E° » E® the results in Fig. Q] suggest differ- 
ent IVR mechanisms for the two states. This can be 
established by correlating the intensities pb a = \(b\a)\ 2 
with the parametric variation of eigenvalues E^fr) i.e, 
the intensity- level velocity correlation function [20. l2l|: 
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where a p and a 



denote the intensity and level velocity 
variances respectively. The parameter r corresponds to 
the resonant coupling strengths in eq. [5] and AE is the 
width of the IVR feature. Recent work 2l| has shown 
that L{,(t) can identify the dominant resonances that 
control the IVR dynamics. In Fig. [5] we show the correla- 
tor Lfc(r) for |1) and |2). Random matrix theory (RMT) 
predicts[20j L{,(t) ~ 0±l/>/N with N being the number 
of eigenstates under the IVR feature and hence ergodicity 
implies a vanishing correlator for any state of choice. It 
is clear from Fig. [2] that several of the correlators violate 
the RMT estimate indicating localization. In particular 
Lfc(r) indicates differing IVR dynamics out of |1) and |2). 
For instance, Lb(k s ff) for the states differ by about 0.2 
which is greater than the fluctuations allowed by RMT 
(« 0.11). Note that the results in Fig.[2]support the local 
RMT approach 2^, 2^1 , developed by Logan, Leitner, and 
Wolynes, which is consistent with the power law decay 
of P b (t) and thus d ~ D (cf. eq.H]). 

We now show that the observed power law in Fig. [T] 
and the slower IVR dynamics of the interior state |2) are 
due to the existence of dynamical traps in the classical 
phase space. First the classical limit Hamiltonian H(I, 6) 
is constructed using the correspondence [2j a,j — ► yije • J 
with (I, 8) being the action-angle variables of Ho. Next, 
classical trajectories with initial conditions such that 
H(1, 6) » E? and actions I restricted to the specific 
state are generated. For every trajectory the dynamical 
function Zk(t) 



y/2I k {t) exp(i6 k (t))_with k 
subjected to the wavelet transform[14j: 



s,f,a,b is 



W g z k {A,B) = A- 1 ' 2 



with A > and real B. The function g(t) — 
(27TCT 2 ) -1 / 2 exp(27riA< — t 2 /a 2 ) is taken to be the sim- 
ple Morlet-Grossman wavelet [14J with A = 1 and <j = 2. 
Eq. [8] yields the frequency content of z k (t) over a time 
window around t = B. In this work we obtain the local 
frequency associated with z k (t) by determining the scale 
(A, inversely proportional to frequency) which maximizes 
the modulus of the wavelet transform i.e., fl k (t = B) = 
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FIG. 2: Intensity-level velocity correlator L^(t) for |6) = |1) 
(circles) and |2) (triangles). The random matrix (RMT) 
estimate is shown by the shaded region. The parameter 
r are taken to be the various resonant strengths and the 
sff, saa, sbb, sfa, sfb, sab, aabb, ffaa, and ffbb resonances 
are denoted by A,B,...,I respectively. The averaging is per- 
formed over a range AE = ±700cm -1 corresponding to the 
width of the IVR feature as seen in the insets (Left inset cor- 
responds to |1) and the right inset to |2)). 



rnaxA\W g z k (A, B)\. This gives the nonlinear frequencies 
n(t) and the dynamics at E = E? is followed in the fre- 
quency ratio spacefill (^s/^a, f2 s /f2b). The frequency 
ratio space is divided into cells and the number of times 
that a cell is visited for all the trajectories gives the den- 
sity plot. We further normalize the highest density to 
one for convenience. Two points should be noted at this 
stage. First, such a density plot is providing a glimpse 
of the energy shell and is reflecting the full dynamics in- 
cluding the important resonances. Thus we are mapping 
out parts of the Arnol'd web i.e., the resonance network 
that is actually utilized by the system. Secondly, we are 
computing a slice of the energy shell and for strongly cou- 
pled systems one expects the phase space structure to be 
different for different slices i.e., nontrivial dependence on 
the angles 6. We thus compute a highly averaged struc- 
ture in the frequency ratio space which is nevertheless 
still capable of providing important information on the 
nature of the classical dynamics. The resulting density 
plots are shown in Fig. [3] for |1) and |2) and look similar 
because E° » E$. 

Fig. [3] clearly shows the heterogeneous or nonuniform 
nature of the density despite angle averaging. This sug- 
gests that at E w E® there are dynamical traps in the 
phase space and hence the dynamics is nonergodic. How- 
ever more important is the nature of these trapping re- 
gions since one expects them to provide insights into the 
IVR dynamics. In Fig. [3] two significant traps corre- 
sponding to il s w 2flf (sff) and another to il a ~ 
(ab) are observed. Note that the f2 a w Of, lock is an in- 
duced effect and in particular persists upon removing the 
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FIG. 3: Dynamical Arnol'd web plotted in the frequency ra- 
tio space (Q a /Q a ,Q a /Qb) for |1) (top) and |2) (bottom). The 
scale for the axes is identical for both the figures. The data 
is obtained by propagating 5000 trajectories with fixed ini- 
tial actions corresponding to the state of interest and varying 
angles such that H(I,6) w E°. The color scale (normalized) 
indicates the number of times that a region is visited by the 
trajectories. Total propagation time T = 10 ps and note that 
the figures look similar already by about T = 2.5 ps. 



geneous nature of the phase space due to resonance trap- 
pings. If the density plots look homogeneous due to the 
absence of any trapping regions then one can associate a 
dimensionality df r = 2 to the frequency space. However 
Fig. [3] show that for both states df r < 2 and one might 
therefore associate a fractal dimension between df r = 1 
and df r — 2. Clearly df r {\) > df r (2) and hence one can 
conjecture that d ~ df r i.e., the effective dimensionality 
of the IVR manifold is the same as the effective dimen- 
sionality of the frequency ratio space or resonance web. 

We conclude by making a few observations. Gambogi 
et al. observed f25J a similar effect in propyne wherein 
the eigenstate-resolved spectra indicated that the combi- 
nation mode v\ + 2vq is much less perturbed by IVR as 
compared to the nearly isoenergetic 3«i overtone state. 
It was argued that such effects are to be expected in large 
molecules. The present example shows that enhanced in- 
stability of overtone states as compared to the combina- 
tion states can occur in few mode systems as well. The 
current study highlights this to be a dynamical effect. 
The decoupling of the (s, /) modes from the (a, b) modes 
for state \2) implies that the full Hamiltonian for E w E® 
is dynamically decoupled into two sub-Hamiltonians: one 
approximately conserving the polyad 2v s + Vf and the 
other conserving the polyad v a +«(,. The precise forms of 
such sub-Hamiltonians is not clear as of now. Finally, the 
extensive afr-lock and the resulting decoupling of the CD- 
bend modes may relate to the observation made by Beil 
et al. on a possible case of an approximate symmetry [13] 
which arises from a near conservation of a formal a' sym- 
metry associated with the bending states. This point, 
however, requires further studies. 



"f a b term from eq. [5] The traps are seen for both states, 
hence the power law behavior of Cb(T) for both states, 
but the extent of trappings is different. The a6-lock is 
more extensive for the state \2) as opposed to the state 
|1). Given the extensive aMock for the dynamics asso- 
ciated with 1 2) one imagines that the CD-bend modes 
get isolated rather quickly from the other two modes. In 
other words, as soon as the energy flows into one of the 
bends the other bend starts to resonantly shuttle this 
energy back and forth resulting in restricted IVR. This 
correlates well with the results in Fig. Q] which shows a 
smaller effective dimensionality of the IVR manifold for 
1 2). Thus one can infer that the restricted IVR for the 
interior state is due to the extensive trapping in the clas- 
sical phase space. The effective dimension of the IVR 
manifold d arising due to a power law behavior of the 
quantum C(,(T) indicated restricted IVR. At the same 
time analysis of the classical dynamics shows the hetero- 
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